Time-invariant person-specific frequency templates in human brain activity 
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The various human brain tasks are performed at different locations and time scales. Yet, we 
discovered the existence of time-invariant (above an essential time scale) partitioning of the brain 
activity into personal state-specific frequency bands. For that, we perform temporal and ensem- 
ble averaging of best wavelet packet bases from multi-electrode EEG recordings. These personal 
frequency-bands provide new templates for quantitative analyses of brain function, e.g., normal vs. 
epileptic activity. 
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Introduction. The various brain tasks (cognitive, 
sensory, motor etc.) are performed simultaneously in 
many locations and operate at different time scales. In 
order to quantify abnormal vs. normal brain function, 
as in the case of epilepsy, it is necessary to develop time- 
invariant templates for characterization of normal behav- 
ior. The challenge is to decipher time-invariant features 
from multi-electrode EEG recordings of brain activity. 
Moreover, from Physics perspective, it is not clear a pri- 
ori that time-invariant elements do exist since the activity 
is inherently nonergodic. 

Here we present a new approach of temporal and en- 
semble averaging of best-bases constructed from Wavelet 
Packets Decomposition (WPD) of the recorded brain ac- 
tivity. The approach is illustrated via the analysis of sub- 
dural EEG (ECoG) recordings from grids of electrodes 
that are placed directly on the surface of the cortex (Fig. 
P). Such recordings are performed to analyze the brain 
activity of epileptic patients under chronic evaluation be- 
fore resection surgery to remove the epileptic focus or foci 
Q Using this approach, we discovered the exis- 
tence of time-invariant personal state-specific frequency 
bands above an essential time scale of about 2.5 min- 
utes. We devised a quantitative measure for comparison 
between WPD bases. Our new person adapted analysis 
can help, for example, in the identification of the epilep- 
tic foci and in the development of quantitative analysis 
methods for early warning of epileptic seizures. As a 
self-consistency test that the frequency bands are not an 
artifact of the analysis, we show that the same templates 
are obtained for subdural and scalp EEG recordings of 
the same person. 

The recorded ECoG signals. The signals analyzed 
here are ECoG recordings from grids of typically ~100 
electrodes. These recordings are obtained from epileptic 
patients undergoing chronic evaluation for surgery. The 
electrodes are spatially distributed over the suspected fo- 
cal region, so that the focus or foci could be localized 
by analyzing the ensemble of signals. The amplitude 
of the signals records the electrical voltage at each elec- 
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FIG. 1: : (a) Subdural EEG (ECoG) grid of electrodes 

placed on the pial surface of the brain for chronic evaluation of 
epileptic patients before surgical resection, (b) Voltage traces 
of ECoG signals. A 10 second time window display of voltage 
traces of 10 electrodes, taken from a multi-electrode recording of 
96 electrodes. 



trode (arguably recording local field potentials [LFP] p|). 
The voltage signals are simultaneously digitized at 112Hz 
(sampling time interval Atmi n — 9 mSec) with a lowpass 
filter up to 40Hz. The analysis is usually performed in 
time windows of Nun — 1024 samples. 
Time- Frequency analysis. In general, the possible 
time intervals for a recorded sequence of Nu n elements 
can range from At m i n = 1 (in units of the sampling time 
interval) to At max = Nu n . In principle, one can ex- 
tract information about N^n time intervals at each of 
the Nun temporal locations along the sequence. How- 
ever, such an N% in matrix for a sequence of only Nun 
data samples must contain redundant information (i.e. 
over-complete representation of the recorded sequence). 
In order to avoid such redundancy, only Nun time- 
frequency locations should be selected, subject to the 
uncertainty constraint between time and frequency reso- 
lutions - At ■ Af = 1. 

Since there are Nun corresponding frequency bands, 
ranging from Af min = 1/N bin to Af max = 1 (in units 
of the Nyquist frequency), each location can be as- 
signed a local relative resolution At/Af out of Nr = 
1 + log2{Nun) possible ratios (for simplicity, Nun of the 
sequences considered here are in factors of 2). It is conve- 
nient to illustrate both constraints as tiling of the time- 
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FIG. 2: An illustration of the Wavelet Packets Decom- 
position (WPD) for a typical signal recorded from a sin- 
gle electrode, (a) A binary tree representation of the best basis 
obtained by WPD. Each node in the tree represents a basis fre- 
quency band. The vertical axis represents the frequency while the 
horizontal axis represents the information cost gain of the basis 
frequency band blocks, (b) The corresponding information dis- 
tribution in the frequency bands of the best basis (darker shades 
represent more information), (c) 2-D tiling representation of the 
WPD of the signal (darker shades represent larger wavelet packet 
coefficients). The vertical axis represents the frequency while the 
horizontal axis represents the time. If the rectangles were colored 
according to their information cost contribution then the informa- 
tion distribution shown in (b) would simply be a sum over the rows 
of the basis coefficients. 



frequency plane with Nun rectangles, each with its own 
aspect ratio (height Af and width At), representing the 
relative resolutions in time and frequency jy, Q • 
The WPD as a binary tree. The Wavelet Packets De- 
composition we use here was devised to partition (tile) 
the time-frequency plane into such rectangles (referred 
to as 'information cells' or 'Heisenberg boxes') pi. Ifj. Ho| . 
Each possible combination of Nun non-overlapping tiles 
that geometrically covers the entire corresponding time- 
frequency plane can serve as a complete basis spanning 
the recorded sequence. 

The WPD is computed by iterating a set of lowpass and 
highpass filters (H and G respectively). The functions 
underlying the expansions of H and G are "wavelets" 
("mother") and "scaling" ("father") functions 11]. At 
each iteration, the wavelet packet coefficients are com- 
puted by convoluting the signal with the filters. Here 
we utilize the WPD using the 'Coiflet' of order 1 
as a "mother" wavelet (smallest time support of all 
'Comets') [II El H El. 

Using the Best Basis algorithm. The WPD gener- 
ates an over-complete representation of the signal. The 
challenge is to select, out of all possible representations, 
the one that is the most efficient in extracting the fea- 
tures of interest. The Best Basis algorithm is a method 
for selecting a basis that spans the signal with a small 
number of significant packets || For that, each 

wavelet packet function is assigned an information cost 
value M q (q) = — q • log 2 (g) where q is the normalized en- 
ergy of the wavelet packet. The total information cost of 
a frequency subband is obtained by summing over all the 
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Viewing the frequency subbands as nodes in a binary 
tree, the selection of the best basis is similar to a binary 
tree search. Starting from the lowest level bands, we se- 
lect for each pair of subbands either the two subbands 
or their joint "parent" band, whichever has the lower 
information cost. The process is repeated at subsequent 
levels, going up the scales, back to the global root. Doing 
so, we select the set of subbands with the lowest possible 
information cost (Fig. |2j). 

Bases Similarity Measure. Every wavelet packets ba- 
sis can be described by the frequency subbands parti- 
tioning and the corresponding information cost of every 
one of the subbands. We suggest using the following in- 
formation cost similarity measure, annotated ICS, for 
comparison between decomposition bases - basisl and 
basis2: 



ICS 



Enl=n2( M nl+M n2 ) 



(2) 



where M n \ and M n2 are the information cost of subbands 
nl and n2 for basisl and basis2, respectively. The sum- 
mation in the nominator is over all the common subbands 
of the two bases. The idea is to compare the information 
cost included in the similar subbands to the total in- 
formation cost. Note, that the measure assumes values 
between (if the bases are totally dissimilar) and 1 (if 
they are exactly the same). 

The Ensemble Best Basis of multi-electrodes 
recordings. In Figs. and |3d we show an example 
of the evaluated Best Bases for the recordings from two 
different electrodes at consecutive time windows (~ 10 
seconds). As can be seen in these figures, the Best Bases 
differ from electrode to electrode (ICS = 0.46 ± 0.20 
between these electrodes) and also vary between con- 
secutive time windows for the same electrode (ICS — 
0.58 ± 0.12 [Fig. Et] and ICS = 0.63 ± 0.25 [Fig. 0b]). 
Looking for invariant Best Bases, we proceed to evaluate 
an ensemble Best Basis for all the electrodes. Following 
previously devised methods 0, 0, 0, ^| , we average 
over the information cost binary tree of each of the L 
recorded signals in the ensemble. This is done by first 
calculating the information cost of all the nodes in the 
binary WPD tree for each of the signals. Next, we eval- 
uate the mean information cost of every node for the L 
signals by simple averaging. Then we apply the Best Ba- 
sis algorithm to the mean values tree. Thus, we obtain a 
basis that may not be optimal for each signal, but rather 
underlines the mean content of the ensemble. The noner- 
godic nature of the brain activity is reflected in the fact 
that the resulting ensemble best basis (EBB) varies be- 
tween successive short time windows (ICS = 0.65±0.33). 
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FIG. 4: The frequency templates - before, during and after 
epileptic seizure, (a) The stability of the frequency template of 
the ECoG recording is demonstrated for time intervals of approx- 
imately 2.5 minutes each, 10 minutes before the seizure onset and 
30 seconds before the seizure onset. During the seizure, the en- 
ergy distribution of the signals changes dramatically |23|. After 
the seizure, the energy distribution changes again, as the energy 
in the high frequencies decreases, (b) Preliminary analysis of scalp 
EEG signals from the same recording of (a) shows the same fre- 
quency template (inter-ictal and post-ictal) and the same temporal 
changes in the energy distribution. 
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FIG. 3: Personal state-specific frequency bands. The fre- 
quency subbands division of the best basis of a single ECoG elec- 
trode signal can change over time, as can be seen in (a) and (b), 
when calculated at 10 consecutive time windows of 9 seconds each 
(we used only 5 levels of decomposition, to avoid a large influ- 
ence due to a negligent number of coefficients, when using short 
time windows). Comparing (a) to (b), it is evident that the fre- 
quency subbands division is also distinct for the different spatial 
locations. Even the ensemble best basis of all the 96 ensemble 
electrodes, as shown in (c), does not yield time- invariant partition- 
ing into frequency subbands. However, temporal coarse- graining 
of the ensemble best basis by averaging over time windows of ap- 
proximately 2.5 minutes produces robust frequency bands that are 
time- invariant for long periods (over 10 minutes). This robustness is 
presented in (d) by examples from the same recordings of such time- 
invariant partitioning into frequency bands (before the seizure) |2l| . 
Nonetheless, the large diversity between individuals, which can be 
observed in the 5 examples in (e), leads us to regard them as state- 
specific frequency bands partitioning or spectral signatures. 



Similarly, if we start with temporal coarse-graining of 
the Best Bases for the individual electrodes, the result- 
ing Best Bases are different from each other and from the 
Ensemble Best Bases. This is a reflection of the inherent 
nonergodicity of the brain activity. 

Time- Invariant Best Bases. However, we did discover 
the existence of an underlying time- invariant (ICS — 1) 



Best Basis in the nonergodic activity. This basis is dis- 
covered by combining ensemble averaging and temporal 
coarse-graining over a new essential time scale. Namely, 
by temporal coarse-graining of the short time EBBs over 
a time window wider than an essential time scale (about 
2.5 minutes) 20]. The latter satisfies the requirement 
that the ICS between EBBs at different time segments is 
larger than 0.95. 

As shown in Fig. the resulting Best Basis is time- 
invariant for recorded periods that are much longer than 
the essential time scale |2lj. Hence, the time-invariant 
bases can be used as a frequency decomposition tem- 
plate for analyzing the recorded brain activity at differ- 
ent times and locations [22^ . 

Personal State-Specific Frequency Bands. Apply- 
ing the new spatio-temporal averaging of the Best Bases 
to the recordings from different persons (we analyzed 
recordings from 12 persons), we found that each has its 
own state-specific time-invariant Best Basis with its own 
characteristic features, as illustrated in Fig. |Hk We em- 
phasize that each of these bases bears resemblance to the 
classical EEG frequency bands, yet has its own specific 
significant deviations from it. Reflecting this notion, we 
note that the inter-patient similarity is ICS = 0.75 ±0.16 
for these examples, which is higher than the similarity be- 
tween the variations of an EBB of a single patient over 
short time windows. However, this inter-patient similar- 
ity is significantly lower than ICS ~ 1, which is measured 
for bases of different time windows of recordings of a sin- 
gle patient above an essential time scale, and could not be 
considered as invariant. We propose that the frequency 
bands are personal specific spectral signatures that can 
be used in patient-specific diagnosis of recorded brain ac- 
tivity. 
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Self- Consistency Test. To further substantiate this 
idea (and that the calculated frequency bands are not an 
artifact), we show a comparison of the frequency bands 
calculated in parallel, for the same person, both from 
ECoG and scalp EEG recordings. As illustrated in Fig. 
EJ between seizure episodes (inter-ictal) the calculated 
frequency bands are almost identical (ICS — 1). We 
also show that the frequency partitioning changes dur- 
ing the seizure (ictal) episodes. Hence, we expect that a 
decomposition of the ictal activity according to the inter- 
ictal bands can help in seizure diagnosis. 
Conclusions. These results illustrate the potential 
value of the personal best-basis to serve as a template 
for quantitative analysis of the epileptic activity 24] . For 
example, the ictal activity can be decomposed according 
to the inter-ictal bands during the chronic monitoring 
of the brain activity and vice versa. Our new person 
adapted analysis can help, for example, in the identifica- 
tion of the epileptic foci. It can also be used to develop 
quantitative person-s pec ific analysis for early warning of 
epileptic seizures |25l l26j. 

Beyond recorded brain activity, we expect the new ap- 
proach to be helpful in revealing the existence of essential 
time scales and time-invariant frequency decomposition 
templates in a wide class of other nonergodic biological 
systems with mult i- time scale dynamics. As in the case 
of the brain, we expect that revealing such hidden tem- 
plates can help in analyzing variations in the systems 
function and performance. 
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